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Abstract. We study the excitation energy for slow changes of the hopping parameter 
in the Falicov-Kimball model with noncquilibrium dynamical mean-field theory. The 
excitation energy vanishes algebraically for long ramp times with an exponent that 
depends on whether the ramp takes place within the metallic phase, within the 
insulating phase, or across the Mott transition line. For ramps within metallic or 
insulating phase the exponents are in agreement with a perturbative analysis for small 
ramps. The perturbative expression quite generally shows that the exponent depends 
explicitly on the spectrum of the system in the initial state and on the smoothness 
of the ramp protocol. This explains the qualitatively different behavior of gapless 
(e.g., metallic) and gapped (e.g., Mott insulating) systems. For gapped systems the 
asymptotic behavior of the excitation energy depends only on the ramp protocol and 
its decay becomes faster for smoother ramps. For gapless systems and sufficiently 
smooth ramps the asymptotics are ramp-independent and depend only on the intrinsic 
spectrum of the system. However, the intrinsic behavior is unobservable if the ramp 
is not smooth enough. This is relevant for ramps to small interaction in the fermionic 
Hubbard model, where the intrinsic cubic fall-off of the excitation energy cannot be 
observed for a linear ramp due to its kinks at the beginning and the end. 



PACS numbers: 71.27.+a, 67.85.-d 
1. Introduction 



In equilibrium thermodynamics, adiabatic processes are denned as quasistatic processes 
without heat exchange with the environment. The entropy remains constant during an 
adiabatic process, while it always increases if the process takes place in a finite time 
and is therefore no longer quasistatic and reversible. These fundamental concepts are 
closely related to the adiabatic theorem of quantum mechanics [1-3] for an isolated 
system which evolves according to the Schrodinger equation 

ihm)) = H{t)\m), (!) 
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with a time- dependent Hamiltonian H(t), i.e., a system that is subject to external fields 
or to changes of its parameters, but not coupled to heat or particle reservoirs. The 
adiabatic theorem states that a system that is initially in the ground state evolves to 
the new ground state during an infinitesimally slow change of the Hamiltonian, whereas 
it cannot follow a parameter change that takes place in a finite time, resulting in a non- 
zero excitation energy. The paradigm for this crossover from adiabatic to nonadiabatic 
behavior in a quantum system is the exactly solvable Landau-Zener model [4,5], i.e., a 
two- level system Hiz{t) = vta z + ^a x that is driven through an avoided level crossing 
with finite speed v > [a z and a x are Pauli matrices). When the system is in the 
ground state |<^o( — oo)) = (1,0) + at time t = — oo, the probability to find the system 
in the excited state \<fti(oo)) = (1,0) + at time t — > oo vanishes exponentially when the 
speed v is small compared to the scale 7 2 /^ set by the gap 7 at the avoided crossing, 



The above Landau-Zener formula can be generalized to various multilevel cases 
[6-9], from which, e.g., the demagnetization probability for the transverse-field Ising 
model was obtained [10]. However, for correlated systems in general the Landau- 
Zener results cannot be directly applied, because essentially all matrix elements of an 
interacting many-particle Hamiltonian change in a complicated way upon variation of 
one of its parameters. The investigation of slow changes of external parameters in 
correlated systems has recently received considerable attention due to its relevance for 
experiments with ultracold atomic gases in optical lattices [11], in which quantum- many 
body systems can be kept under well-controlled conditions. In those systems time- 
dependent control of the parameters is not only of practical importance (as discussed 
below), but it also allows to test fundamental theoretical predictions. For example, the 
Landau-Zener result was indeed experimentally confirmed in a Bose-Einstein condensate 
loaded into an accelerated optical lattice [12]. 

Various slow parameter changes in many-body systems have recently been studied 
[13-31]. For a general ramp the system is initially in the ground state and some 
parameter of the Hamiltonian is then changed to a new value within a time interval r, 
either linearly or nonlinearly with time. To investigate the crossover from the extreme 
nonadiabatic limit r = (i.e., a sudden quench of the Hamiltonian) to possibly adiabatic 
behavior in the limit r — > 00 a measure for the degree of nonadiabaticity is needed. A 
popular quantity for this purpose is the excitation energy AE(t) after the ramp, i.e., 



where Eq(t) is the ground-state energy of the Hamiltonian after the ramp. For initial 
states at non-zero temperature, the entropy increase provides a more natural measure 
of nonadiabaticity in general. However, entropy is uniquely defined only for thermal 
equilibrium, and thus it can only be computed after the ramp is complete and the system 
has thermalized. On the other hand, isolated many-body systems do not necessarily 
thermalize quickly after changes in the Hamiltonian [32-45], in particular for integrable 
systems, as demonstrated experimentally with ultracold gases [46]. In contrast to the 



\(i/)(t -> oo)||0i(oo)))| 2 ~ exp(-7r 7 2 /^). 




(2) 
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entropy the internal energy is always well-defined, regardless of whether the system 
passes through a series of thermal or nonthermal states in the limit of a quasistationary 
process. 

In the present work we only consider systems that are initially in the ground state. 
If the excitation energy AE(t) vanishes in the limit of long ramp times, r — > oo, the 
system is considered to behave adiabatic. It is expected that the excitation energy is 
still small for finite ramp times r, just as in the Landau-Zener formula, when the ground 
state is protected by a gap for all parameters throughout the ramp [24]. However, the 
excitation energy is not exponentially small (AE(t) oc exp (const /r)) in general. As we 
will show below, the asymptotic decrease of AE(r) for large r can depend both on the 
intrinsic properties of the many-body system and on the ramp protocol. In particular 
for gapped systems the ramp protocol can be used to make AE(t) arbitrarily small, 
but in general it often vanishes only algebraically. This behavior is known from the 
Landau-Zener model, where the excitation is exponentially small only when the avoided 
level crossing is traversed from t = — oo to t = +00, whereas the excitation probability 
is proportional to 1 /r 2 and hence much larger if the evolution takes place from t = to 
t — 00, i.e., starting exactly at the center of the level crossing [18,47]. 

The situation is completely different for gapless systems, such as the exactly solvable 
one-dimensional transverse-field Ising model in which the gap vanishes at exactly one 
value of the transverse field. When the magnetic field is ramped across this critical 
point the excitation energy is [14-16] 



with a rational exponent r] = ^. Similar results were obtained for a number of other 
quantum critical systems, such as the Bose-Hubbard model [21] or the random field 
Ising model [20]. However, the existence of a quantum critical point is not a necessary 
condition to obtain a nonanalytic relation AE(t) [23,24]. Equation Q, with various 
values of the exponent 77, holds for ramps within gapless phases of several gapless 
systems [24]. For a continuous bath of harmonic oscillators, which model the low- 
energy excitations of a large class of systems, the exponent rj for a slow squeeze of the 
oscillator mass depends on the spatial dimension [24]: An analytic relation AE(t) ~ 
t~ 2 is found for all dimensions d > 3, while r\ is noninteger for d = 2. For d — 1, the 
thermodynamic limit does not commute with the limit of large r, i.e., the prefactor in 
Eq. ([3]) increases with system size [24], suggesting that adiabatic behavior is impossible 
for that class of one-dimensional systems. 

The excitation energy during a nonadiabatic ramp, and its dependence on the 
ramp duration r is not only a fundamental property of a quantum many-body system, 
but it is also of practical interest for experiments with cold atomic gases. Various 
ramping procedures are used in experiment to transform one phase into another, and 
the available time for the process cannot be too long in order to avoid extrinsic losses. 
On the other hand, whether theoretical predictions are actually observable in experiment 
can depend in a subtle way on the unavoidable excitation during the preparation of the 




(3) 
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state [17,48-53]. When the ramp duration is fixed to a given maximum value, it thus 
becomes important to find the optimal ramp through which a given point in parameter 
space can be reached through minimal excitation of the system [25]. In general, it is 
plausible that any additional term in the Hamiltonian should be switched on slowly, so 
as to build up the correlations that it favors without incurring to high energy cost, and 
increasing the speed at later times. 

In view of these issues the question arises to what extent the dependence of the 
excitation energy on the ramp duration r is determined by intrinsic properties of the 
system, and to what extent it is influenced by the details of the ramp. In this paper 
we give a perturbative argument that holds in the limit of small ramp amplitudes and 
allows to separate an intrinsic contribution to the excitation energy and a ramp shape 
dependent contribution. In some cases the latter can mask the intrinsic contribution 
such that the behavior of the excitation energy in the limit of long ramp times r is 
completely determined by the ramp shape. Furthermore, we present results for the 
excitation of the Falicov-Kimball model after various ramps. In this model, which can 
be solved exactly using nonequilibrium dynamical mean-field theory (DMFT), Eq. ([3]) is 
found to hold with an exponent r\ that is different for ramps across the metal-insulator 
transition, within the metallic phase, and within the insulating phase. Our numerical 
results for rj in this model support the scenario obtained from the perturbative argument. 

The paper is organized as follows. In Sec. [2] we show results for the excitation 
energy in the Falicov-Kimball model in nonequilibrium dynamical mean-field theory. In 
Sec. |3] we develop a perturbative argument for small ramp amplitudes and discuss the 
implications for gapped and gapless systems, such as the metallic and Mott insulating 
phases of the Falicov-Kimball and the fermionic Hubbard model. A conclusion in Sec. |4] 
closes the presentation. 

2. Ramps in the Falicov-Kimball model 



Below we present results for the excitation energy AE(r) for ramps of different types 
in the Falicov-Kimball model [54], with Hamiltonian 



Here q and f> are annihilation (creation) operators for the itinerant and immobile 
electrons, respectively, and n\ = cjcj (n{ = ) are their local densities. Hopping 
between sites i and j, with amplitude Vij(t) = V(t)tij, is possible only for the mobile 
c particles. Note that although the / electrons are immobile, the equilibrium state of 
-^Falicov-Kimball does not correspond to one quenched / configuration but rather to a state 



2.1. Model 





(4) 
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with annealed disorder, where each / state contributes according to the free energy of 
the c particles. 

In the context of dynamical mean- field theory (DMFT) [55], which becomes exact in 
infinite dimensions [56] , the Falicov-Kimball model has a long history because it can be 
mapped onto a solvable single-site problem [57-60]. A Mott metal-insulator transition 
occurs at a critical interaction U c for half-filling (at density n c = rif = |), as well as a 
transition to a charge-ordered state at sufficiently low temperatures. The physics of the 
Falicov-Kimball model thus partly resembles that of its parent, the fermionic Hubbard 
model, 

^Hubbard (t) = ^ ^i^ylaCja + U^^n^Un - fi ^ ( 5 ) 

«?>=T,| * i,c=T,J. 

with two mobile spin species. 

DMFT can be applied to nonequilibrium situations [39,45,61-68], in which case 
the effective single-site problem for the Falicov-Kimball model is still quadratic and can 
be solved using equations of motion [62]. Here we extend the exact solution of the 
Falicov-Kimball model for an interaction quench [39] to a numerical solution that can 



be applied to arbitrary time dependencies in V(t) and U(t) ( Appendix A ). This allows 
us to study the excitation after ramps of the hopping or the interaction strength. We 
employ a set of hopping amplitudes for which the density of states has a semielliptic 
shape, 

p( e ) = ^E 5 ( e -^) = ^^ 3 ^' ( 6 ) 

k 

where are the eigenvalues of the hopping matrix tij and L is the number of lattice 
sites. Furthermore, we consider only the homogeneous phase at half-filling, for which 
the chemical potential is fixed at fi = U/2 and the /-orbital energy at Ef = 0. In this 
case the critical interaction for the equilibrium Mott transition is U c = 2V [58]. 



2.2. Linear ramp protocol 

We consider linear ramps in the Falicov-Kimball model (Sj) in DMFT for the 
homogeneous paramagnetic phase at half-filling. We assume that the system is in the 
ground state for times t < 0. For < t < t the hopping parameter V is changed 
according to the ramp protocol 

f Vi t<0 

V(t) = < Vi + AV r(t/r) 0<t<r (7) 
V f = Vi + AV t>r, 

where Vi is the initial hopping amplitude, r is the total ramp time, AV is the ramp 
amplitude, and r(x) is the ramp shape. The latter is a monotonously increasing function 
with r(0) = and r(l) = 1. We set the energy scale by V, = V = 1, so that time is 
measured in units of 1/V. (From now on we set K — 1.) The energy of the system per 
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Figure 1. Excitation energy ([2]) after linear ramps of the hopping parameter [Eq. 0, 
r(x) — x] within the metallic phase (U — 1, Vf — 2), within the insulating phase 
(U = 3, Vf = 0.5), and across the metal-insulator transition (U = 1, Vf =0). The 
energy scale is set by V% = V = 1. The curves become independent of r in the quench 
regime r < l/V. The solid black lines, with a slope 1/2, 1, and 2 (from top to bottom), 
correspond to the asymptotic behavior ©. Inset: Internal energy E(t) [Eq. ([5])] during 
ramps (0 of the hopping amplitude in the Falicov-Kimball model (r(x) = x, U = 1, 
and Vf = 0), using various ramp durations r. For t < and t > r, the energy is 
constant. The solid black line is the internal energy E(t) in the ground state at U = 
1 and hopping V(t). 



lattice site is given by 



v(t) ■£ ty(t)(ci(t) Cl (t)) + u(t) J2(4(thWHt)f, (*)> 



The excitation after a ramp is then obtained from the difference 02]), where Eq(t) is the 
energy ([8]) of the ground state of the Hamiltonian after the ramp. The DMFT solution 



for ramps in V(t) (and also U(t)) is described in Appendix A 



The time evolution of the energy ([8]) is plotted in the inset of Fig. [T] during a 
ramp (J7)) with linear profile r(x) = x. For small ramp-durations (r = 1), the energy 
rises linear with time. In this case, the system is essentially quenched, i.e., its state 
\i}(t)) remains unchanged during the ramp, and the energy is thus only determined by 
the ramp protocol, E{t) w (ip(0)\H(t)\tp(0)). In the opposite limit r — > oo, the energy 
adiabatically follows the ground-state energy E (t) for hopping parameter V(t) (solid 
line in inset of Fig. [1]), in accordance with the adiabatic theorem. 
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We now focus on the excitation AE{r) after the ramp, which is plotted in Fig. [T] 
for linear ramps ([7j) within the gapless metallic phase (U = 1, Vf = 2), within the 
gapped insulating phase (U = 3, Vf = 0.5), and across the metal-insulator transition 
(U = 1, Vf = 0), which occurs in the equilibrium system at U = 1 and V = 0.5. 
From Fig. [T]one can estimate the crossover timescale r quen ch, which separates the regime 
in which the state of the system cannot follow the parameter change (r < r qucnch ) 
from the adiabatic regime in which AE(t) decreases with increasing ramp-duration r 
(t > ^quench)- Independent of the ramp parameters, r qucnc h turns out to be of the order 
of few times the inverse bandwidth. The decrease of AE{r) for r > r quenc h can be fitted 
with a power law (EJ) for r > 10. The exponent turns out to be a rational number, which 
depends only on the phase in which the system is before the ramp (metallic phase for 
U < 2, insulating phase for U > 2) and after the ramp (metallic phase for U < 2Vf, 
insulating phase for U > 2Vf). These results can be summarized as 

{t~2 linear ramp across the transition, 
t -1 linear ramp in metallic phase, (9) 
r~ 2 linear ramp in insulating phase. 

How do these exponents arise and how do they depend on the ramp shape? 
Further data show that the exponent i] = \ for the excitation across the metal-insulator 
transition is independent of the ramp shape r(x). At present we have now simple 
explanation of this exponent. It would be interesting to determine how this exponent is 
related to the critical behavior of equilibrium correlation functions, such as the density of 
states at the transition [59] . On the other hand, the behavior for ramps within either the 
metallic or the insulating phase will be explained in the next section by a perturbative 
argument, which applies to small ramps of arbitrary shape in any quantum system. In 
particular we will see that the exponent rj = 1 is a consequence of the non-Fermi-liquid 
behavior of the metallic phase in the Falicov-Kimball model, while the exponent t] = 2 
in the insulating phase is not an intrinsic property of the Falicov-Kimball model but is 
in fact due to the linear ramp shape. 

3. Small ramps of arbitrary shape without traversing phase boundaries 

Our numerical results for ramps of the hopping amplitude in the Falicov-Kimball model 
show that the exponent rj in Eq. (Q does not depend on the precise values of the ramp 
parameters V, arid Vf, but only on the thermodynamic phase of the initial and final state. 
This finding suggests to study the excitation energy perturbatively in the limit of small 
ramp amplitudes, but for arbitrary ramp shapes and ramp durations. In the remainder 
of this section we will derive the excitation energy AE{r) up to second order in the 
ramp amplitude for an arbitrary Hamiltonian. In particular we discuss the asymptotic 
behavior of AE(t) in the limit r — > oo and how it may be influenced by the ramp shape, 
and illustrate these general results with data for the specific case of the Falicov-Kimball 
model. 
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3.1. Perturbative result for the excitation energy 

We consider the general Hamiltonian 

H{t) = H + K (t)W, (10) 

where H is the Hamiltonian before the ramp, W is the operator that is switched on, 
and n(t) is the ramp function. As in Eq. (JTj), we characterize n(t) by the ramp amplitude 
An, the ramp duration r, and the ramp shape r(x), i.e., n{t) = An r(t/r). In order 
to expand AE(r) for fixed ramp duration r and ramp shape r(x) in powers of An, 
we decompose the quantum state \ip(t)) of the system in the instantaneous eigenbasis 
|0n(£)) °f the Hamiltonian (fTOl) . which satisfies the condition 

H(t)\<f> n (t)) = €n(t)\<f>n(t)) (11) 

at any instance of time. We assume the \4> n (t)) to be nondegenerate for simplicity. 
After fixing the phase of the eigenvectors in a convenient way we obtain the eigenstate 
decomposition of \ip(t)) as 

\m) = Y, a ^ elf ° dsen(s) \^> ( 12 ) 

n 

so that the Schrodinger equation implies 

ijfnit) = £ e^* 6 ™^ (M)\j t \K{t)), (13) 

m 

using the notation e nm (t) = e n (t) — e m (t). The matrix element on the right-hand side 
of Eq. ffl3l) is given by 

€nm(t) (Mt)\j t \Mt)) = (Mt)\H(t)j t - ^-H{t)\(j> m {t)) 

= (0n(t)|^|0 m (t)) 

= AK r -^^(<f> n (t)\W\(f> m (t)), (14) 
r 

where the first equality follows from Eq. (fTTI) and the last from the explicit form of the 
Hamiltonian [Eq. ffTOD ]. 

Because the system is assumed to be in the ground state |0o(O)) of Hq for t < 
0, Eq. ( Ti~3l) must be solved with the initial condition a m (0) = <5 m o- Together Eq. ( fl4l) 
this implies that a n (t) = O(Ak) for n ^ 0. In order to obtain the leading term in the 
expansion of a n (t) (for n ^ 0) we can thus restrict the sum in Eq. (fT3l) to the single 
term m = 0, 

..w = a- f«-iMMW»t,^«(.) + 0(Ak2) (1B) 

Jo r e n0 (t) 

This expression was used before as a starting point for the discussion of ramps across a 
quantum critical point [14]. Here we study ramps which do not cross a phase boundary, 
and we assume that the instantaneous eigenenergies e n0 (t) and eigenfunctions \(f> n (t)), 
which depend on time t only through the parameter k, can be expanded around k = 0. 



Near-adiabatic parameter changes in correlated systems 



9 



Since a n (t) is already of order O(Ak), e n0 (t) and |</> n (^)) in Eq. ([15]) can be replaced 
by e„ = e„ (0) and \(j) n ) = \cj) n (0)), respectively The excitation energy, AE{t) = 
T En^o^oW \an(t)\ 2 , is then given by 

AE(t) = An 2 £{t) + 0(Ak 3 ) (16) 

= f°° — R(u)F(ut) (17) 
Jo v 



£(r) 



R{") = y^|(0n|W|0o)|^(^ 



L 



F(x) 



ds r "^" lxs 



s e° 



(19) 



Eqs. (I16I)- (I19I) constitute the main result of this section. The correlation function 
R(uj), which can be interpreted as the spectral density of possible excitations induced 
by the operator W, is independent of the ramp shape r(x) and the ramp duration r. 
Conversely, the ramp spectrum F(x) does not depend on the Hamiltonian but only on 
details of the ramp. For continuous ramp shapes r(x) it follows that F(x) — > for |x| 
— > ±oo, such that F(ut) becomes increasingly peaked around u = in the limit r — > 
oo. In fact, making the replacement F(ut) oc 5{uj)/t is equivalent to Fermi's Golden 
Rule for |a n (t)| 2 , and the nonadiabatic excitation (|17|) is due to deviations of F{ujt) 
from 5{uj). The crossover scale r quonc h that was discussed in Sec. [2] is thus given by the 
value of r below which F(ut) = F(0) + 0{uj 2 t 2 ) is approximately constant over the 
entire bandwidth Q of R(u), i.e., r quench pa 1/Q. In Sec. 13.51 we confirm this estimate 
numerically for small interaction ramps in the metallic phase of the Falicov-Kimball 
model. 

In the following we will analyze the asymptotic behavior of Eq. (TTTI) in the adiabatic 
limit, r — > oo. For this we need the behavior of the ramp spectrum F(x) at large values 
of x, which follows from Eq. ffl9|) as 



™/(x) ; w . th /(x) ^oo 

x a 

where the exponent is given by a = In if the nth derivative of r[x) is discontinuous 
(i.e., the (n — l)st derivative has a kink), but all lower derivatives are continuous; this 
behavior follows from the Riemann-Lebesgue lemma [69] . For example, in case of a linear 
r(x) = x the first derivative r\x) = Q(x)Q(l—x) is discontinuous at x = and x — 1, so 
that the ramp spectrum F(x) decays like x~ 2 [cf. Eq. (129 &l) below]. In general, when the 
ramp shape has a finite number of kinks, the large-x asymptotics of the ramp spectrum 
([19]) is given by a finite sum of oscillating terms, F{x) ~ | J2k fk cos(ukx) + 5k\ 2 /x a . By 
choosing a smooth ramp one can always increase the exponent a or even make F(x) 
decay exponentially for x — » oo. However, in practice ramp protocols often have kinks 
that lead to a power- law decay (|20|) . 

To estimate the magnitude of the integral (TTTI) in the limit r — ► oo we distinguish 
two cases, namely (i) the gapless case, in which R(u) vanishes like a power law at 
u — 0, and (ii), the case of a gapped excitation spectrum, in which R(uj) has a finite 
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gap r2g a p above uo = 0. In both cases we assume that R(u) is zero beyond some high- 
frequency scale Q, although the argument remains valid if R(u>) vanishes exponentially 
for ijj > Q. Since J °° du R(uj) = -^[(0 o |W /2 |0o) — (^ol^l^o) 2 ] we can also assume that 
any singularities of R{oj) are integrable. 

3.2. Case (i): Gapless excitation spectrum 

In this paragraph we discuss the case in which the excitation spectrum R{uj) is gapless 
and vanishes like a power law at uo = 0, 

R(u) ^°uj u with v > 0. (21) 

First we assume a > u, where a is the exponent that characterizes the ramp shape 
[Eq. (l2~0"i) ]. Writing R(cu) = u> u R(ou), the integral ( fT71) becomes, after a change of 
variables, 

£{t) = \ [ di/- 1 E(ir 1 )F(x). (22) 
rU Jo 

Using the asympotic behavior (120]) we find that the integral in this expression remains 
finite in the limit r — > oo, so that in this case the leading contribution to the excitation 
energy is given by 

a>u: E(r) T ^^- = E intI (r), (23) 

with C = R(0) J °° dxx u ~ 1 F(x). It is important to note that the exponent does not 
depend on the ramp shape, but only on the density of possible excitations above u — 0. 
Because the latter is an intrinsic property of the system we will refer to ^i I1 tr(' r ) a s the 
intrinsic contribution to the excitation energy in the following. In principle, a ramp 
between two parameter values can always be made so smooth that E- in tr{^) becomes the 
dominating contribution to the excitation energy (i.e., a > u), as in Eq. (1231) . However, 
as we will see in the following paragraph, if the ramp is not smooth enough (i.e., if 
a < v), the intrinsic contribution will be masked by a nonuniversal contribution that is 
essentially determined by the ramp shape. 

For the case of a < v we estimate the integral (JTTj) as follows. For the moment we 
assume that the spectral density R(u) has no singularities at finite frequencies and use 
the bounds 

u u d -u)< R(u) < u u C 2 6(fi 2 - u), (24) 

with positive constants C\, C 2 , fli, fl 2 - Together with Eqs. (l20ll and (JTTj) we obtain for 
the excitation energy for r — > oo, 

a < v. ^ < £(t) < — (25) 

log(rfl) log(r^) 
a = z/: — — — < £{r) < C 2 — — — , (26) 

with positive constants C[ and C 2 . The upper bound holds because f(x) = 0(1) in 
Eq. (1201) . To obtain the lower bound it is sufficient to note that although f{x) can 
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have infinitely many zeros, the moving average f{x) = f* + dxf(x) over any small 
finite interval of given length h is larger than some positive constant. This property is 
satisfied in particular when the ramp shape has a finite number of kinks, as discussed 
below Eq. (1201) . Finally we note that Eqs. (123j) and (1261) hold also if R(u) has integrable 
singularities, because a small frequency interval around each of them contributes to 
the integral (ITTj) in the same way as the gapped spectrum [Sec. 13. 3j . namely oc r~ a 
[Eq. flZZD]. 

The result that is stated in Eqs. (1231) and (126"]) has a simple interpretation: Kinks 
in the ramp shape increase the probability of excitations to high energy states, as 
expressed by the slowly decaying tail of the ramp spectrum F(x). When the ramp 
is not smooth enough the integral (ITTj) is therefore dominated by the high-frequency 
part of R(uj), leading to a nonuniversal, ramp-shape dependent excitation energy. For 
the often- considered linear ramp (a = 2) any intrinsic contribution Si ntr (T) with v > 2 
will therefore be unobservable in S(r). This is precisely what happens for weak-coupling 
interaction ramps in the Hubbard model, as discussed below in Sec. 13.51 

3.3. Case (ii): Gapped excitation spectrum 

We now turn to the case of an excitation density which has a gap f2 gap at uo = 0. The 
integral (ITT)) then starts at the finite lower bound fi gap , such that F(x) can be replaced 
by its asymptotic behavior ( 1201) in the entire integration range. As a consequence we 
have 

£{T) ™ JL f" dx gMZM . (27) 

a 'gap 

The integral gives a finite constant in the limit r — > oo provided that R(u) is not 
singular. Otherwise the integral may give a r-dependent but bounded contribution, as 
shown in the next subsection for ramps in the insulating phase of the Falicov-Kimball 
model. The gapped case [Eq. (127]) ] is thus similar to the gapless case with a < v 
[Eq. (1251) ]. In both cases the excitation energy is dominated the high-frequency behavior 
of F(x) and is therefore completely determined by the ramp shape, while the intrinsic 
contribution ( [23]) is unobservable. 

Our analysis so far can be summarized as follows. The excitation energy after 
a ramp may be either dominated by the intrinsic contribution (123]) or set by ramp- 
shape dependent terms [Eqs. ( 125]) . ( 1261) . and ( 1271) ]. depending on the large-frequency 
asymptotics ([20]) of the ramp spectrum and the small-frequency behavior (121]) of the 
excitation density. This fact will be illustrated in the following two subsections for ramps 
in the insulating and metallic phase of the Falicov-Kimball model and the Hubbard 
model. 

3.4- Insulating phase 

In previous subsection we have shown that the excitation energy after a ramp within 
a gapped phase behaves in a nonuniversal way because the intrinsic contribution ( [23] ) 
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Figure 2. Ramps shapes given by Eqs. (|28a[) - (|28cj) (left panel), and corresponding 
ramp spectra [Eqs. (|29a[) - (|29c[) ] (right panel). The Fresnel oscillations in F n (cu) are 
due to the discontinuity in the derivatives of r n (x) at x — and x = 1. 



vanishes. A significant dependence of the nonadiabatic excitation energy on the ramp 
shape is therefore expected also for ramps with finite amplitude. In the following we 
will demonstrate this fact for ramps within the insulating phase of the Falicov-Kimball 
model, where it turns out that the asymptotic behavior for r — > oo is indeed correctly 
described by the analytic expression ( !27|) that was obtained for small ramps. 
For this purpose we focus on three particular ramp shapes, 

r\{x) = x (28a) 
1 — cos(7rx) 



r 2 (x) 
r 3 (x) 



TTX 



COS [TlX) Sill TTX I 



(286) 



:28 c) 



7T 



COS (CV) 



Here r n (x) is chosen in such a way that its nth derivative is discontinuous at x = and 
x — 1 (Fig. [2^), i.e., r' n (x) oc sin n (7nc) for < x < 1. The corresponding ramp spectra 
[Eq. (USD] are 

F 2 (x) 
F 3 (x) 



7r 4 1 + cos(o;) 2 



2 

32tt 4 



7T- 



U 2 ) 2 

- cosfo;) 



(29a) 
(296) 

These functions vanish like F n (x) ~ x~ 2n for x — > oo (Fig. [2b). We now perform ramps 
of the hopping amplitude V{t) = V i J r (Vj — Vi)r n (t/r) in the Falicov-Kimball model, 
with Vi = V = 1 as the energy scale. We consider only the paramagnetic insulating 
phase at half-filling, i.e., U > 2 = 2V* and U > 2Vf. 

The excitation energy after such ramps is plotted as a function of the ramp duration 
r in Fig. [31 The curves can be fitted with power laws (131) for large r, with an exponent 
7] = 2 for the linear ramp (128 aft and rj = 4 for the ramp (12861) . respectively. For the 
ramp (128 cl) the results are consistent with an exponent t] = 6, but the excitation energy 
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Figure 3. Same as Fig. [TJ for ramps within the insulating phase. r(x) = r n {x) 
[Eqs. (|28a|) - (|28cl) ]. U = 3, Vf = 0.5. The black linear lines correspond to power-law 
behavior ([3]) with r\ = 2 (n = 1) and t] = 4 (n = 2), and rj = 6 (n = 3). 



is too small for a power-law fit in the accessible range. Hence the large-r behavior of 
the nonadiabatic excitation in the case of ramps with finite amplitude turns out to be 
the same as in the limit of small ramps, i.e., a power law with an exponent that is 
determined by the singularities of the derivatives of the ramp shape [cf. Eq. ff27j) ] rather 
than by intrinsic properties of the system. 

To check Eq. ( IPTI) explicitly for arbitrary ramps we would have to compute the 
density of excitations R{ui), for which no general solution is available. Nevertheless one 
can derive an expression in the atomic limit and compare the resulting excitation energy 
to ramps deep in the insulating phase. For ramps of the hopping V(t), the operator W 
in Eq. ( Fl8|) is given by the kinetic energy operator. In the case of half-filling for both 
mobile and immobile particles there is exactly one particle per site in the ground state 
for V = 0. Therefore each hopping process creates exactly one doubly-occupied site, 
and the function R(u) consists of a single delta peak at u = U. When R(u) oc 5{uj — U) 
is inserted into Eq. fflTl) one obtains 

AE(t) oc F(Ut), (30) 

i.e., the Fresnel oscillations in ramp spectra such as Eqs. (129 aft to (129 aft become visible 
in the dependence of the excitation energy on r. This result should only be slightly 
modified for ramps deep in the insulating phase (U ^> V) , assuming that the delta-peak 
is then only slightly broadened and shifted in position. In fact, as seen in Fig. H] the tail 
oscillations of Fi(x) [Eq. (129 ap . Fig. [2] are also apparent in the excitation energy for 
ramps between states with U ^> V. In the limit r — > oo these oscillations are washed 
out because R(u>) has a finite bandwidth AQ for V > 0, such that the integral (fT7|) 
averages over many oscillations for r ^> 1/AQ. 




0.1 1 10 



I 

Figure 4. Excitation energy [2] after the linear ramp (|28a[) of the hopping amplitude 
within the insulating phase [U = 10, Vi = 1, Vf = 0.5], compared to Eq. (|50"|) for small 
ramp amplitudes [F± given by Eq. (|29aj) ]. 



3.5. Metallic phase 

As an application of Eq. ( ITTjl to ramps in a gapless phase we consider the turn-on of 
the interaction in the Falicov-Kimball model and the Hubbard model. For the following 
discussion it is convenient to write the Hamiltonian in momentum space 

H = ]>> fc(T - ^)ci a c ka + U(t)D (31) 

ha 

D = J2 H ^ = C l+q| C fc| C L-qT C fc'f ( 32 ) 

i k,k',q 

Furthermore we change the notation with respect to Eq. (j3J) to allow for a unified 
description of the Hubbard model, where both spin species are mobile (efcj = e*.;, /if = 
/if ) and the Falicov-Kimball model, where we take spin f to be immobile (e^ = 0, /if = 
fi — Ef, /i| = /i). We consider ramps at half filling = U/2) in which the interaction 
is changed from zero to a finite value, U(t) = AUr(t/r). 

For a ramp of the interaction strength in the Hubbard model and the Falicov- 
Kimball model, the operator W in Eq. (TT0T) is given by the double occupation (|32l) . As 



shown in Appendix B , the excitation density R(u) at U = can be expressed in terms 
of the second-order contribution to the self-energy 

R(u) = -- Im££We qi + z0). (33) 



7T 

— a;<e q | <0 



For comparison to our DMFT results we evaluate Eq. fl33l) in the limit of infinite 
dimensions [56] where the self-energy is independent of momentum q, [70] and the 
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qr-summation can be replaced by an integral over the density of states P|(e), 

1 f° 

R(u)) = / dep t (e) Im sf } (e + u + iO) . 




(34) 



For the Hubbard model, the second-order self-energy is given by [70] 




(u)= / dfip(fi-iu) / dup(u)p(ji-u) ~ |p(0) 3 cu 2 . (35) 



Jo Jo 



In accordance with Fermi liquid theory the imaginary part of the self-energy vanishes 
oc u 2 , thus leading to well-defined quasiparticle excitations in the metallic phase of 
the Hubbard model. On the other hand, the imaginary part of the mobile electron 
self-energy in the Falicov-Kimball model remains finite at uj = due to the scattering 
off fixed impurities. Its value can be obtained easily from the exact solution of the 
Falicov-Kimball model in DMFT [60], 



Here nj is the average density of localized particles, such that n/ = 0.5 in case of half- 
filling. Eqs. (l3~5l) and ( l36l) can then be inserted in Eq. ( l34l) . which in turn determined 
the intrinsic component (|23|) of the excitation energy, 



As discussed above, the intrinsic contribution can be masked completely by a ramp- 
shape dependent contribution if the ramp is not smooth enough, i.e., when the exponent 
in Eq. (I2"0"j) satisfies a < 1 or a < 3 in case of the Falicov-Kimball and Hubbard model, 
respectively. However, the discussion below Eq. fl20l) shows that F(x) decays at least 
oc x~ 2 if the ramp is continuous, i.e., if it does not contain any abrupt finite changes. 
Hence we conclude that the intrinsic component ( 1381) is always dominant for ramps that 
turn on a small interaction in the Falicov-Kimball model. This is consistent with our 
numerical results for ramps in the metallic phase of the Falicov-Kimball model which 
are not shown here, namely, that the 1/t behavior for the metallic phase in (Q does 
not only hold for linear ramps (Fig. [1]), but for all three ramps (I28aj) to ( 128 cj) . 

The situation is very different for ramps in the Hubbard model. Because the 
intrinsic contribution (!37|) vanishes oc r~ 3 for r — > oo it is negligible with respect to 
the high-frequency contribution ( 1251) for linear ramps, where the ramp spectrum decays 
as x~ 2 [Eq. (1 2 9 <a[) ] . This is consistent with results of Mockel and Kehrein [31], who 
computed the excitation energy after a linear ramp of the interaction in the Hubbard 
model using Keldysh perturbation theory and found AE(r) ~ t~ 2 for r — > oo. 

The exact expression for R(u) at all frequencies in the Falicov-Kimball model at 
U = [Eqs. ( 134"1) and (1361) ] allows us to evaluate S(t) for arbitrary interaction ramps 
U(t) = AUr(tr) at finite ramp durations, and compare to the DMFT result for AE(t). 
In Fig.[5]this comparison is done for a smooth Gaussian ramp [inset (a) in Fig. [5], which 



Im £y. M H = (1 - n f )n f p(uj) ^ (1 - n f )n f p(0). 



(36) 




(37) 
(38) 
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Figure 5. Excitation AE(t) in metallic phase of the Falicov-Kimball model, after 
Gaussian ramps [Eq. c 2 = 36, see inset (a)] from interaction U = to AU. 

The curve AU = is obtained from Eq. (|17|) . where R(u>) is obtained from Eqs. (|3^|) 
and (|36[) . and -F(u;) is given by Eq. ([59]) . The excitation spectrum (solid line) is 

plotted in the inset (b), together with 0.03 F(lut) (dashed lines) to illustrate the effect 
of r (see text). 



we define by 

r'{x) = c x exp [ - c 2 (x - 1/2) 2 ], (39) 

F(x) = —± exp - — for c 2 > 1, (40) 
c 2 V 4c 2 / 

where ci is a normalization constant to satisfy the condition J^dxr'(x) = r(l) — r(0) 
= 1, and c 2 is chosen such that r'(x) is sufficiently small at the boundary x = and 
x = 1, i.e., the expression for F(x) holds up to terms which are exponentially small in c 2 . 
For U < 1 numerical results for the excitation energy (scaled with the ramp amplitude 
AU 2 ) agree very well with the analytical expression ( TT7T) . evaluated using Eqs. ( 1341) . 
(]36p . and PU|) . This corroborates the validity argument of Sec. [3] and shows that it 
provides correct estimates for the nonadiabatic excitation energy after ramps which are 
not too large in amplitude. 

Inset (b) of Fig. [5] illustrates the origin of the crossover from small to large 
ramp times r. For fast ramps, e.g., r = 1, the ramp spectrum F(ujt) averages over 
the entire bandwidth Q of the excitation spectrum R{oS). As a consequence, the 
excitation energy becomes independent for ramp times smaller than the quench time 
scale r < 1/Q = r qucnch . Indeed we see in the numerical data that the weak-coupling 
quench time scale fits with the estimate 1/Q ~ 1/4. For larger quench times, e.g., 
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t > 10 in inset (b) of Fig. \5[ the ramp spectrum F(ur) probes only the linear small-u; 
behavior of R(w), leading to the universal power-law behavior AE(t) ~ 1/r for r > 10. 

4. Conclusion 

We presented a general perturbative analysis of the excitation energy due to slow ramps 
of a parameter in a quantum system without crossing of phase boundaries, motivated 
by our numerical results for the Falicov-Kimball model obtained with nonequilibrium 
dynamical mean-field theory. We demonstrated that the excitation energy vanishes 
algebraically for large ramp duration r [Eq. ([3])] under rather general circumstances. 
The exponent r) can depend on one hand on the spectrum of the correlation function of 
the operator that is switched on, and on the other hand on the differentiability of the 
ramp function. Which of these influences dominates in rj depends on the low-energy 
behavior of the excitation spectrum compared to the spectrum of the ramp protocol. In 
practice, any experimental ramp protocol can always be considered as differentiable on 
a short enough timescale. Our conditions on the degree of differentiability have to be 
interpreted in the sense that a ramp protocol must be considered as not differentiable 
if the slope or any higher derivative changes on a timescale shorter than the inverse 
bandwidth of the system. 

For ramps in gapped systems the asymptotic behavior of the excitation energy 
depends only on the ramp protocol and can be made as small as desired by use of 
increasingly smooth ramp shapes. By contrast, for ramps in gapless systems the low- 
energy excitation spectrum has no effect on rj if the ramp is not smooth enough. Only 
if the ramp is sufficiently smooth does r] become ramp-independent and reflects the 
low-energy excitation spectrum of the system. For the fermionic Hubbard model this 
implies that a linear ramp from U = to a small value of U leads to an unnecessarily 
large excitation energy with rj — 2, which can be reduced to the intrinsic exponent rj = 
3 if the ramp shape has at least two continuous derivatives. 

Our results also indicate that in the Falicov-Kimball model the exact expression 
for the excitation energy in the limit of small ramps provides a good estimate up to 
quite large ramp amplitudes. This suggests to use the perturbative expression, which 
is valid for arbitrary systems, as a guide for finding ramp protocols that connect fixed 
parameters of the Hamiltonian and minimize the excitation energy for a given ramp 
time, thereby improving the preparation of states in experiments with ultracold atomic 
gases. 
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Appendix A. Solution of the Falicov-Kimball model in nonequilibrium 
using DMFlfl 

In this appendix we describe in some detail how the Falicov-Kimball model (j4]) with 
arbitrary time-dependent hopping amplitude V(t) or time- dependent interaction U(t) 
is solved using nonequilibrium DMFT. In DMFT, local correlation functions of the 
lattice model are obtained from an effective impurity model in which a single site 
is coupled to a self-consistently determined environment [55]. The mapping of the 
lattice model onto the single-site model, which becomes exact in the limit of infinite 
dimensions [56], can be formulated either in imaginary time, which yields a theory 
for thermal equilibrium, or in real time (using the Keldysh technique), which yields a 
theory that can be applied to various nonequilibrium situations [61]. For the Falicov- 
Kimball model the action of the single-site model can be reduced to a quadratic one [57], 
such that nonequilibrium correlations function can be determined from a closed set of 
equations of motion [39,62-67]. 

Because we are interested in the transient time evolution of a system which is in 
thermal equilibrium for times t < (i.e., its initial state is given by the density matrix 
p oc exp[— f}H(0)]), we use contour-ordered Green functions with time arguments on the 
contour C that runs from t — to t max (the maximum simulated time) on the real axis, 
back to t = 0, and finally to t = —if} along the imaginary axis [71]. The local Green 
function is then given by 

G(t,t') = -i^Tr[e-^T c c(t)c\t% (A.l) 

where Tc is the contour ordering operator, and c^(t) are annihilation (creation) 
operators of the mobile particles in the Heisenberg picture with respect to the time- 
dependent Hamiltonian. Up to a factor i the imaginary-time Green function of the 
interacting equilibrium state is recovered from Eq. (1A.1I) when both time arguments 
are on the imaginary-time portion of the contour. On the other hand, when both time 
arguments are on the real branch we obtain the real-time Green functions G < (t,t') 
= iTr[e-P H( -V$(t')c(t)} and G > (t,t') = -iTr[e-^ H ^c(t)c^(t% from which various 
thermodynamic observables can be calculated. In particular, the internal energy per 
lattice site (1A.2I) is given by 

E(t)=d t G < (t,t')\ t=t ,, (A.2) 

which follows directly from the equations of motion of the lattice system, assuming 
spatial homogeneity. 

The DMFT equations for the Falicov-Kimball model with time-dependent 
interaction were stated in detail in Ref. [39] and the appendix of Ref. [72], where an 
analytical solution for the case of a sudden switch of the interaction parameter is given. 

| In the appendices, r denotes imaginary time, not to be confused with the ramp duration r in the 
main text. 
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The local Green function flA.lh is determined from the equations of motion [39], 

G(t, t') = w Q(t, t') + Wl R(t, t'), (A.3) 

[id t + fi]Q(t, t') - [A * Q] (t, t') = 6 c (t, 0, (A.4) 

[id t + U(t)]R(t, - [A * R] (t, t') = 5 c (t, 0, (A.5) 



where W\ = 1— Wq denotes the average number of localized particles (which is fixed in the 
homogeneous phase), and A(t,t') is the coupling to the environment which is obtained 
by a self-consistency condition. The product [A * B] (t, t') denotes the convolution of two 
functions along the contour C, and 5c(t,t') is the contour delta function. Throughout 
this paper we use the semielliptic density of states (jSJ), in which case the self-consistency 
condition takes the simple form [39] 

A{t,t') = V 2 G{t,t'). (A.6) 

Equations (1A.3I) to (1A.6I) form a closed set of integro-differential equation on C. 
In this paper we consider also the case of a time-dependent hopping amplitude in (j4]), 
i.e., we assume Vij(t) = V(t)tij, where the hopping amplitude V = V(0) sets the 
energy scale and the ty are dimensionless. This case can be mapped onto the case 
of an Hamiltonian with time-independent hopping and time-dependent interaction in 
the following way: The action exp[— % J dtH(t)} of the lattice model is invariant under 
a simultaneous scaling of the Hamiltonian H(t) = VH(t)/V(t) and transformation to 
new time variable t(t) = f dt'V(t') /V. By definition, H has time-dependent interaction 
U(t) = UV/V(t) but constant hopping V, such that Eq. flA.6jl is valid. Under the 
same transformation of time variables, the hybridization function A(t,t') transforms as 
A(ti,t 2 ) = V(ti)A(ti, t2)V{t<2) /V 2 . Hence the change of the time variable leads to a 
replacement of the self-consistency equation ( 1A.6I) by 

A(t,t') =V{t)G{t,t')V{t'). (A.7) 

By discretizing of the contour C, Eqs. (IA.3I) . flA.41) . (IA.5I) . and (IA.7I) can in principle 
be reduced to the inversion of a matrix whose dimension is given by the number of mesh 
points along C [62,63]. This is approach is not suitable here, because it would require 
a infinite length of the contour in case of initial states at zero temperature (j3 — > oo). 
A different approach first parametrizes contour Green functions in terms of various real 
and imaginary time components, and uses Langreth rules to derive separate integral 
equations for each component [71]. One can then remove the imaginary time branch by 
a partial Fourier transform to Matsubara frequencies and analytical continuation to real 
frequencies [68]. The solution of integro-differential equations such as (1A.4I) and (1A.5I) 
in this way is described in detail in Ref. [68]. In the following we will therefore only 
briefly restate these equations to mention the differences that arise from the fact that 
we are not solving a single equation (such as Eq. ( 1A.4I) for given A), but a nonlinear set 
of equations. 

In the following we adopt the notation of Ref. [68] . When both time arguments of 
a contour Green function A(t,t') are on the imaginary time portion of the contour, we 
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obtain the Matsubara component which can be represented in the form 

A (-vr, _ iT ') « £ e ^'-^a(icu n ). (A.8) 

n 

The function a(u) can be continued to real frequencies. Using this parametrization to 
rewrite Eqs. (IA.4I) and (IA.5I) [68] , we obtain the well known cubic equations for the local 
Green function of the homogeneous phase in the Falicov-Kimball model [59] 

#(o; + zu) w + /i _y(o) 2 ^(a; + zO) a; + -T/(0) 2 ^(o; + i0)' 1 j 

The solution of this equation (with negative imaginary part) describes the initial state. 

Next we consider the retarded component A r (t,t') = Q(t - - A < (t,t')] 

of the contour Green function, for which Eqs. (1A.3I) . flA.4j> . (1A.5j) . and (1A.7[) read [45] 

G r {t, t') = w Q r {t, f) + w 1 R r {t, t'), (A.10) 
[td t + fi}Q T {t,t') - J dsV{t)G T {t,s)V(s)Q T (s,t') = 0, (A.ll) 

[id t + fi-U]R T (t,t') - f dsV(t)G T (t,s)V(s)R T (s 1 t > ) = 0. (A.12) 

Jt' 

These equations must be solved for t > t' with the initial condition G r (t,t) = R r (t,t) = 
Q T (t,t) = —i. In contrast to the case discussed in Ref. [68], this is a set of nonlinear 
integro-differential equations. However, the equations are causal in the sense that the 
differential is always determined by an integral over the function at earlier times. Hence 
the solution of Eqs. flA~T0|) to flATT]) is very similar to the solution of an ordinary 
differential equation, and the presence of a nonlinearity does not lead to additional 
difficulties. 

In addition to equations for the retarded and Matsubara components one has to 
consider equations for the mixed components A~^(t, r) = A(t,—ir). In terms of the 
partial Fourier transform 

A^{t,iu n )= drA^{t,r)e^ nT , (A.13) 



we obtain, after analytically continuing iui n — > u± = u ± iO, 

GT{t,u)±) =w Q"(t,cu ± ) + w 1 R~"(t,uj±), (A.14) 

[id t + {j]Q^{t,u±)- [ dsV(t)G v (t,s)V(s)Q^(s,uj ± ) 
Jo 

= V(t)G^(t,u ± )V(0)q M (u ± ), (A.15) 

[id t + n -U}R^ (t, u±) - [ dsV(t)G I (t,s)V(s)R^(s,uj ± ) 

Jo 

= V(t)G^(t,u ± )V(0)r M (u ± ), (A.16) 

to be solved with the initial condition A _1 (0, u±) = ia M (u±). Finally, the lesser 
component satisfies 

G < (t,t')=w Q < (t,t')+w l R < (t,t'), (A.17) 
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[^ + /i]g < (t,t')- / dsV{t)G I {t,s)V{s)Q < (s,t') 



Jo 




(A.18) 



[2ft + /i (t, t') - / dsV(t)G I (t,s)V(s)R < (s,t') 



Jo 




(A.19) 



to be solved with the initial condition A < (0, t) = A~(0, t) for A = Q, R. Together with 
the hermitian symmetry, A a (t,t') = A r (t',t)* and A~^(t,u) = A~(u*,t)* for Q, R, and 
G the set of equations is closed. 

Note that the expansion of Eqs. (TOP . fTA~4j) . (TO) , and flA~7j) into Eqs. (IA~T0D - 
( 1A.19I) has the following technical advantage. The numerical solution of Volterra integro- 
differential equations can be easily implemented such that the time discretization error 
scales as e ~ At p = [t/N) p for At — > 0, with p > 1 [73]. To study the excitation 
energy ([3]) for large r from the difference ([2]), the absolute energy (1A.2[) . must be 
determined with an relative accuracy of the order of t~ v . Assuming the above error 
scaling, sufficient accuracy of the energy with respect to AE(t) thus requires N ~ 
t ~i+v/p timesteps for r — ► oo. On the other hand the Green function must be stored 
at 0(N 2 ) time points and hence N is the limiting numerical factor (we go up to A" w 
10000). It is thus crucial to use an algorithm which is correct up to high-order in At 
when the exponent t] is large. For ramps in the insulating phase, e.g., rj > 2 is found 
[cf. Eq. (J9J)], and these results could not be obtained using the lowest-order trapezoid 
approximation (p = 1), but we used higher-order schemes instead (p = 5). 

Appendix B. Excitation density in the noninteracting limit of the Hubbard 
and Falicov-Kimball model 

For a ramp of the interaction strength in the Hubbard model and the Falicov-Kimball 
model, the operator in Eq. fTTOT) is given by the double occupation ( 132|) . and the excitation 
density R{u) can be evaluated at U — 0. For this purpose we introduce the imaginary- 
time ordered correlation function 



where the expectation value (■) = Tr[e _/3H -]/Tr[e~^ ff() ] is taken in the noninteracting 
state at temperature T = 1/(3 (the limit T — ► is taken at the end), and D is given 
by (!32l . (T T is the imaginary time ordering operator, and D{r) are Heisenberg operators 
with respect to Hq). Because D contains an even number of Fermi operators, the 



R(t-t') = (T t D(t)D(t>)}) 



(B.l) 
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Figure Al. Left panel: Diagrammatic representation of Eq. (|B.3|) . Lines represent 
the noninteracting momentum-resolved Green function g^ a (iuj m ) = l/(iuj m — e qrT ) for 
a =t (solid lines) and a =[ (dashed lines). Momentum is conserved at the vertices, 
frequency irj n enters at the left vertex. Right panel: Transformation of the Matsubara 
sum (|B.5[) to the real frequency interval (|B.6|) . using the usual expression Yliu ^(^m) 
= (i/3/2ir) § c dzf(z)w(z), where f(z) is the Fermi function and w(z) is some analytic 
integrand. The integrand in Eq. (|B.5l) has a branch cut at z = —ir\ n due to the branch 
cut of along the real axis, and a pole at e q -\. Then the contour G\ is transformed 
into C2, which yields (|B.5|1 . using that the Fermi function is periodic under shift with 
bosonic Matsubara frequencies, f{u> — irj n ) = f{ui). 



function R(t) satisfies periodic boundary conditions on the imaginary time contour 
r G [0, j3] and can be expanded in bosonic Matsubara frequencies r\ n = 2nn/T, R(ir) n ) = 
J Q dTR(T)e~ irinT . Using the Lehmann representation one can show that the excitation 
density ( fT8l) may be obtained from the unique analytical continuation of R(ir) n ) from 
rj n > to the upper half of the complex frequency plane, 

R(cu) = 1mR{u + i0). (B.2) 

7T 

To calculate R, the expectation value (1B.1I) is factorized using Wick's theorem and 
transformed to bosonic Matsubara frequencies. It turns out that the only nonvanishing 
contractions for r] n ^ is given by 

fci,fc2,q r,s,m 

x 9l 2 i(iu s )g k ^(iuj r )gl 2 _ q _ ki l (iuj m + iu s - iu r + irj n ), (B.3) 

where g° (iu m ) = l/(iu m — e qa ) is the noninteracting Green function at momentum q, 
and iuj m are fermionic Matsubara frequencies. The expression has a simple diagrammatic 
representation (Fig. lAlb ). The diagram is split into one Green function line g^(iu m ) and 

(2) 

the remainder, which we identify as the second-order contribution S J to the selfenergy, 
R(iVn) = Y^gq^iuimjE^^irjn + iu m ). (B.4) 

q iw m 

One can now transform the Matsubara summation into a frequency integral, where it 
must be taken into account that the self-energy has a branch cut along the real 
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axis with S(c<j ± iO) = =)=Im E(u;) (Fig. lAlb ). The result is 

(B.5) 



,(2),. , V \ A fl J mS qTV + i0 ) 



UJ -lTj n - e qT 



9 

where /(e) is the Fermi function. Eq. (IB .51) is finally continued to the real frequencies 
by replacing ir\ n — > uj + iO, and the spectrum (IB. 2ft is obtained as 

R ( u ) = I + " /(^i)] Im S 9l ( w + ( B - 6 ) 

9 

Taking the limit of zero initial temperature then yields Eq. (1331) . 
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